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ABSTRACT 


Using  a  series  of  19  explosions  with  accurately  known 
epicenters  within  a  500  km2  area  of  the  Nevada  Test  Site,  the 
location  effectiveness  is  demonstrated  of  applying  predetermined 
travel-time  anomalies  to  a  limited  network  of  teleseismio.  sta¬ 
tions  (comprised  of  between  4  and  13  stations  greater  than  1900 
km  distance).  Three  different  travel-time  tables  were  used: 
Jeffreys-Bullen;  Herrin  1961  version;  and  Herrin,  November  1966 
version;  and  two  different  computer  programs:  LOCATE  and  SHIFT, 
the  former  which  minimizes  the  sum  of  squares  of  residuals  and  ’ 
the  latter  which  minimizes  the  sum  of  squares  of  relative  resi¬ 
duals.  The  mean  location  error  for  the  19  known  epicenters,  ob¬ 
tained  without  time  anomalies,  is  about  26  km,  and  with  anomalies 

as  less  than  3  km,  regardless  of  travel-time  table  and  regard¬ 
less  of  program. 

It  is  further  demonstrated  that  neither  the  number  of  sta¬ 
tions  in  the  range  of  3  or  4  to  13  nor  the  distance  aperture  of 
the  network  has  an  effect  on  the  location  of  known  surface  events 
although  the  azimuth  aperture  does. 


Confidence  estimates  are  made  in  three  ways:  the  standard 
confidence  ellipses;  maximum-relative-error  polygons;  and  standard- 
deviation  contours  about  the  final  solution.  It  is  shown  that  by 
applying  travel-time  anomalies,  the  standard  confidence  ellipses 
which  estimate  the  reliability  of  the  data  in  a  least  squares  ’ 
sense,  can  be  reduced  in  area  by  factors  of  1/5  to  1/152  and 
still  enclose  the  true  epicenter. 

.  A  discU5si°n  given  of  the  s*nbility  of  travel-time  anom¬ 
alies  across  the  Nevada  Test  Site  area,  and  of  some  problems  in¬ 
volved  in  determining  usable  anomalies  from  earthquakes. 
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INTRODUCTION 


It  has  been  well-established  that  station  time  residuals 
and  travel-time  anomalies  (relative  residuals)  are  not  constant 
but  change  significantly  from  one  epicentral  region  to  another 
at  a  rate  dependent  upon  the  region  and  the  station  or,  with 
travel- time  anomalies,  upon  the  station-pair  separation  (see, 
for  example,  Chiburis  and  Dean,  1965;  Chiburis,  1966a, b,  and 
1968).  However,  the  time  errors  are  reasonably  well-behaved 
within  each  region  such  that  they  can  be  predicted  (time- 
calibrated)  for  additional  events  occurring  in  the  same  region. 

Any  location  scheme  which  assumes  that  there  are  no  time 
errors  other  than  those  due  to  reading;  or  that  a  regional  or 
azimuthal  correction  is  valid  for  all  stations  in  a  network; 
or  that  a  single-station  correction  for  one  region  is  adequate 
for  another  will  necessarily  do  a  poor  job  of  locating  epicen¬ 
ters.  This  report  then,  is  primarily  concerned  with  the  tele- 
seismic  location  accuracy  obtainable  for  a  particular  region 
by  using  either  station  residuals,  travel-time  anomalies,  or 
no  corrections  at  all.  Secondarily,  we  are  concerned  with  the 
performance  of  two  computer  programs  presently  in  use  at  the 
Seismic  Data  Laboratory,  LOCATE  and  SHIFT.  The  principal  dif¬ 
ference  between  the  two  programs  is  that  SHIFT  minimizes  in  a 
least-squares  sense  the  relative-anomaly  errors  rather  than 
the  absolute  residuals. 

Also  investigated  is  a  technique  whereby  limits  of  the 
solution  are  estimated  either  on  the  basis  of  an  acceptable 
network  standard  deviation  of  time  errors  or  of  maximum 
relative- time  errors  at  any  station  pair  within  the  network. 
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DEFINITIONS 


Travel-Time  Anomaly. 

The  travel-time  anomaly  at  station  i  relative  to  station  j 
is  defined  as 

Ai/j  Ti  -  Tj  -  Hi  +  H.  +  e.  -  e. 

where  T  is  the  observed  arrival  time,  H  is  the  expected  travel 
time  according  to  some  travel-time  and  distance  relationship, 
and  e  is  a  correction  for  ellipticity.  In  this  report,  values 
of  H  were  computed  from  the  Jeffreys-Bullen,  Herrin  1961,  and 
Herrin  1966  (November)  travel-time  tables.  The  above  definition 
of  anomaly  makes  no  assumptions  concerning  crustal  P-wave  veloc¬ 
ities  (for  station  elevation  corrections);  it  simply  measures 
the  net  effect  of  all  causes  of  travel-time  errors  at  any  sta¬ 
tion  with  respect  to  any  arbitrary  reference  station. 

Residuals. 

The  residual  at  station  i  is 

Ri  =  Ti  -  To  ‘  (Hi  +  V 

where  Tq  is  event  origin  time  and  d±  is  a  correction  both  for 
ellipticity  and  for  the  elevation  of  the  ith  station  assuming  a 
value  for  the  angle  of  incidence  and  for  the  P-wave  velocity 
between  the  station  and  sea-level.  Here,  H  is  computed  from 
the  Herrin  1961  table  which  is  the  standard  relationship  in  the 
particular  version  of  LOCATE  used  for  this  report. 
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METHOD 


The  program  LOCATE  uses  the  standard  Geiger  technique, 
briefly  described  as  follows: 

Let  the  errors  which  are  to  be  minimized  in  a  least-squares 
sense  be  defined  as 


where  R?  is  the  observed  residual  at  station  i,  R.  is  the 
previously-determined  residual  for  that  region  and  R!  is  a  theo¬ 
retical  residual  such  that  1 


R 


f 

i 


3  H 

Z>V  =  air 


3H . 

dX  +  jt  d(p  + 


3H.  3H. 

rr  dz  +  -?r  dT0 


where  A,*,  Z,  and  Tq  are  respectively  the  event  longitude,  lat¬ 
itude,  depth,  and  origin  time,  and  H  is  the  expected  travel  time. 

As  a  change  in  origin  time  is  the  same  as  a  change  in  expected 
travel  time, 


Also,  to  compare  the  two  programs  more  closely,  depths  were  re 
strained  throughout  such  that 


Therefore 

3H  3H. 

Ri  =  8 X—  M  +  dT0 

T  .  N 

t  is  desired  that  be  a  minimum  for  N  stations  so 

3  N  2 

djr  (  .1  Ei>  1  0  k  =  1,2,3 

K  1=1 

where  p.  =  dA,  =  d$,  and  y,  =  dT  . 

o  o 
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This  differentiation  yields  three  normal  equations  which  can  be 

solved  simultaneously  for  the  errors  in  the  event  parameters  dA ,  dd>, 

and  dT  . 
o 

The  program  SHIFT,  on  the  other  hand,  defines  an  anomaly  error 
at  station  i  relative  to  station  j  as 

dAi/j  =  Ai/j  -  h/i 

where  is  the  previously-determined  relative  anomaly  for  that 

region.  The  relative  errors  to  be  minimized  in  a  least-squares 
sense  are 

Ei/j  =  dAi/j  ‘  dAi/j 

where  dA^.  are  the  theoretical  relative  anomalies  such  that 

dAi/j  ±  ^  «  +  d*  +  <=j 

where  e.  is  the  average  error,  or  bias,  at  the  reference  station  j. 

It  is  desired  to  minimize  J  E?  .  ,  so,  as  before, 

i=l  1/3 

^  Ei/j)  1  0  *  *  1-2.3 

where  Ml  =  dA ,  u2  =  d*.  and  U3  =  c, .  Solving  these  three  normal 
equations  yields  the  event-parameter  corrections  dA ,  d*,  and  the 
reference-station  bias  c.. 

1 
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DESCRIPTION  OF  THE  DATA 

All  of  the  time  data  used  in  this  study  were  derived  from 
nuclear  explosions  detonated  within  the  Nevada  Test  Site  (NTS) 
area  (Figure  1).  Table  I  lists  the  event  information.  The 
events  selected  as  references  from  which  residuals  and  anomalies 
were  measured  are  Bilby  or  Tan,  the  series  of  Bronze,  Corduroy, 
and  Buff,  or  the  series  of  Nash,  Agile,  and  Commodore.  These 
several  events  were  necessary  so  that  residuals  or  anomalies 
could  be  obtained  for  an  adequate  number  of  recording  stations. 
Table  II  lists  the  station  anomalies,  relative  to  RK-ON,  com¬ 
puted  from  the  three  travel-time  tables,  in  addition  to  the 
residuals  computed  from  the  Herrin/ 19 61  table.  A  key  letter, 
indicating  which  series  of  events  were  used  as  references  in 
determining  the  corrections,  is  given  in  the  last  column. 

PROCEDURE 

Station  records  of  all  explosions  were  routinely  read,  with 
the  identical  networks  of  stations  and  arrival  times  being  used 
for  both  programs.  Depths  for  all  events  were  restrained  to  the 
surface.  Raw  arrival  times  (Table  III)  were  input  to  SHIFT  and 
used  within  the  program  in  conjunction  with  the  appropriate  in¬ 
put  travel-time  anomalies.  For  LOCATE,  station  arrival  times 
were  corrected  by  the  residuals  prior  to  their  input  into  the 
program.  These  input  times  for  LOCATE  are  obtained  by  sub¬ 
tracting  the  residuals  in  Table  II  from  the  arrival  times  input 
to  SHIFT  (Table  III). 

A  description  on  the  use  of  program  SHIFT  is  given  in 
Appendix  I. 


RESULTS  OF  LOCATION 


Without  Time  Corrections. 

neithTable  Z  HStS  the  l0Cati°n  Srr0rS’  in  when 

times^thrt  alV°r  an°malieS  are  applied  t0  event  arrival 
ZZn  h  ^  8  be"  C°rrected>  ■  for  station  ele- 

LOCATE  1PtlClty,‘  The  ">eon  error  for  17  events  using 

LOCATE  is  seen  to  be  25.8  to  and  using  SHIFT  (Herrin  61)  25  9  tan 

the  results  fro,  the  two  programs  using  the  sate  travel-tim^ 

tables3"6558"'  ^  aSraamant-  Wth  «>•  -d  Herrin  66 

’  8  r8SUlts  usin8  SHIrT  are  27.7  tan  and  20.8  tan  re- 

apeotrvely.  It  is  not  Known  at  this  time  if  the  apparently 

e  er  results  obtained  by  using  the  Herrin  66  table  are  sig- 
m  leant  or  not.  It  would  be  necessary  to  locate  a  larger 
sa;ple  of  events  to  determine  the  effectiveness  of  this  table 
e  ne  works  used  in  this  study  and  for  the  NTS  area. 

the  be"':  reSUltan"  aTOrage  Srrors  of  about  26  tan  demonstrate 

re!.d  h°Pe  ^  d°  "han  n°  aiLowances  are  made  for 

uals  (LOCATE)  or  anomalies  (SHIFT). 

With  Time  Corrections. 

When  the  residuals  or  anomalies  are  determined  for  a 

Itltion  h-T-0"’  SUCh  ^  "eVada  TeSt  Site’  and  f°r  each 
atation  which  is  to  be  used  in  subsequent  location  networks 

a  resultant  location  errors  can  be  reduced  by  at  least  an 
order  of  magnitude.  Table  V  lists  the  location  errors,  in 
l  ometers ,  when  either  residuals  (LOCATE)  or  relative  ano¬ 
malies  (SHIFT)  are  applied  to  the  event  arrival  times.  The 
mean  error  for  17  events  using  LOCATE  is  now  seen  to  be  2  98  km 

fromTtnh  SHIFT  <Herrin  61>  2-8G  Agai">  tha  values 
from  either  progrem  are  in  agreement,  although  individual  event 

This  diff  1  ar-by  33  mUCh  3S  4-1  k”  between  the  two  methods. 

outing  enCef-”P  ’  °f  C0Ur5e>  th”  the  P"««.  are  com¬ 

puting  in  significantly  different  ways  but  yield  about  the 
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same  answer  on  the  average.  Using  the  J-B  and  Herrin  66  tables, 
the  SHIFT  mean  errors  are  2.92  km  and  2.59  km  respectively  which 
suggests  that  when  anomalies  are  applied,  any  reasonable  travel- 
time  table  is  adequate. 

Actual  Time  Errors. 

Although  the  events  used  for  locating  in  this  study  were  of 
a  reasonable  size  which  made  film  reading  straightforward,  time 
errors  at  the  stations  due  to  misreading,  imperfectly-known 
anomalies,  etc.,  were  far  from  negligible.  Actual  time  errors  at 
station  i  relative  to  mean  zero  for  the  jth  event  are  shown  in 
Table  VI,  computed  as 


0  =  a3  _  t 

/ o  i/r  i/r 


where  A  is  the  anomaly  for  the  jth  event,  A  is  the  previously- 
measured  input  anomaly  for  the  region  of  the  jth  event,  and  E^ 
is  the  mean  error  for  the  jth  event.  Also  included  in  this  table 
are  the  errors  for  A  =  0  (which  assumes  that  there  are  no  ano¬ 
malies).  Event  standard  deviations  are  given  at  the  far  right 
of  the  table.  As  shown  in  the  table,  the  actual  standard  devi¬ 
ations  when  anomalies  are  not  used  average  about  0.86  sec,  with 
errors  as  high  as  1.9  sec.  When  allowances  are  made  for  the 
anomalies,  the  actual  standard  deviations  average  about  0.16  sec, 
a  reasonable  figure  for  reading  error  alone,  although  individual 
errors  are  as  high  as  0.8  sec.  The  point  to  be  made  is  that  the 
set  of  events  used  for  testing  the  validity  of  applying  anomalies 
is  not  unusual,  because  reading  or  other  time  errors  are  not 
particularly  small,  some  even  being  quite  large. 


A  similar  study  is  presently  being  undertaken  in  which  the 
set  is  composed  of  low-yield  events,  at  least  in  a  signal-to- 
noise  ratio  sense.  However,  selecting  this  set  so  that  the  events 
are  recorded  at  a  significantly  smaller  size  than  those  used  in 
this  study  is  difficult  because  several  of  the  events  listed  in 
Table  I  have  quite  low  recorded  signal  amplitudes  due  to 
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deliberately-reduced  station  magnifications.  The  computed 
magnitudes,  where  known,  are  given  in  the  following  list  for  all 
17  events: 


Event 

Magnitude 

Fore 

5.2 

Buff 

5.1 

Chartreuse 

5.  ? 

Auk 

4.9 

Piledriver 

5.5 

Bourbon 

5.1 

Dumont 

5.5 

Agile 

Not  Calculated 

Nash 

5.2 

Commodore 

5.7 

Greeley 

6.3 

Klickitat 

5.0 

Turf 

4.9 

Piranha 

Not  Calculated 

Scotch 

5.5 

Corduroy 

5.6 

Bronze 

5.2 

Film  Reading  of  First  Extrema. 

One  of  the  methods  used  in  this  study,  which  can  substan¬ 
tially  reduce  errors,  is  to  permit  readings  of  times  other  than 
those  for  first  motion.  Regardless  of  event  size,  first  motion 
is  usually  difficult  to  read  consistently  to  within  0.1  sec. 
However,  by  reading  the.  arrival  time  of  the  first  extremum 
(peak)  rather  than  first  motion  at  all  stations,  timing  is  con¬ 
siderably  more  precise.  Therefore,  most  of  the  events  in  this 
study  were  read  for  this  phase,  although  a  few  were  so  large 
the  traces  went  off  scale  forcing  the  reading  of  first  motion. 
The  excellent  location  results  show  that,  by  reading  either 


first  extrema  or  first  motions,  one  does  not  lose,  and  may 
actually  gain,  network  capability. 

.vumber  of  Stations. 

The  number  of  stations  of  which  the  networks  were  comprised 
for  locating  the  set  of  17  NTS  events  varied  between  4  and  13. 
Figure  2  shows  the  location  errors,  in  kilometers,  as  a  function 
of  the  number  of  stations  for  SHIFT-H61,  both  with  and  without 
anomalies.  The  results  without  anomalies  show  no  dependence  on 
the  number  of  stations,  and  those  with  anomalies  are  too  few  to 
derive  a  clear  relation.  This  result  suggests  that  it  is  not 
necessary  to  require  a  large  number  of  stations  for  locating 
accurately  if  proper  allowances  are  made  for  the  station  travel¬ 
time  anomalies. 

Effect  of  Aperture. 

Figure  3  shows  the  location  errors  in  kilometers,  of  the 

17  events  as  a  function  of  azimuth  aperture  which  is  calculated 
as 

de.  =  |e  -  e  .  I 

i  1  max  mm 1 

where  0  is  epicenter-to-station  azimuth.  Included  on  the  figure 
are  the  results  obtained  both  with  and  without  travel-time 
anomalies.  The  errors  appear  to  have  a  fairly  certain  dependence 
on  the  aperture  of  the  azimuth,  especially  for  the  results  ob¬ 
tained  with  anomalies.  The  actual  dependence  is  hard  to  ascertain 
due  to  lack  of  data  points  at  apertures  less  than  80°.  This 
effect  is  presently  being  investigated. 

The  effect  of  distance  aperture  on  the  location  error  is 
negligible,  as  seen  in  Figure  4,  both  with  and  without  anomalies. 


CONFIDENCE  REGIONS 

mannerTfidenr?-reSi°nS  "ere  C°mpUted  within  SHIFT  in  the  usual 

trav  1  t6'8'’  lnn’  1965>  ^  311  events’  both  Mith  and  without 
travel-time  anomalies.  Table  VII  lists,  by  event,  the  areas  of 

the  computed  ellipses  and  the  factors  of  reduction  in  ellipse 

when  anomalies  are  used.  Without  anomalies,  two  events 

level  With  anomalies,  all  events  were  within  the  ellipse  and 
ellipses  were  reduced  in  area  by  factors  of  5  to  152,  with 
an  average  reduction  of  about  45.  This  reduction  in  ellipse 
area  points  out  the  necessity  and  value  of  travel-time  anomalies. 

Maximum  Relative  Errors. 

.  “hen  computati°ns  are  made  for  confidence  regions,  the 
variances  or  standard  deviations,  of  the  final  solutions  are 
used  in  the  estimates.  These  variances  are  merely  indicators 
of  the  goodness-of-fit  in  the  least-squares  procedures  and  if 
the  number  of  degrees  of  freedom  is  small,  the  variances  and, 
ce,  the  confidence-region  estimates  become  unrealistic, 
n  aot.when  as  few  as  three  stations  are  used,  no  estimate 
at  all  can  be  made  of  the  location  error,  as  the  solution  is 
unique).  For  example,  the  events  Bourbon  and  Scotch  have 

(TlblCtVTely’a°tUal  Standard  Nations  of  0.08  and  0.23  sec 
(Table  VI,  With  travel-time  anomalies)  and  N=4  and  5.  The 

1  9*“  T?  T  BOUrb°n  “  7118  ta”2  a"d  thS  1U°ation  “«*  ^ 

1.9  km,  while  for  Scotch  the  area  is  2309  km2  and  the  error  is 
^  ^•11B°th,ellipSe  ab"°™ally  high  due  strictly  to 

data  fr“  T"  "  time 

as  thl  d  t  !  thS  "atW°rk  apa«^s  are  as  good 

the  data  and  apertures  from  the  other  networks.  Therefore, 

when  few  stations  are  used  in  locating  events,  it  may  be  of  value 

to  estimate  the  lotion  error  by  considering  the  maximum  permis- 

r-tJ  th  h1”6  errOT  “  aly  Stati°n-Pair  in  the  network 

61  a"  3  statl5tical  estimate  using  the  variance  of  the 
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goodness-of-fit.  In  this  way  fairly  good  estimates  of  location 
error  can  be  made  even  when  the  number  of  stations  is  three. 

Figure  5  shows  the  maximum  relative  errors  (X100)  using 
Bourbon  data  for  a  grid  of  locations  about  the  final  SHIFT 
solution.  This  grid  is  computed  as  an  option  within  SHIFT  by 
subroutine  SIGRID.  The  scale  in  latitude  is  1.0  km  and  in 
longitude  2.4  km  between  center  points  of  each  number  field  of 
4  spaces.  With  this  proportion,  the  scale  is  equal  in  all 
directions.  The  errors  at  each  grid  point  are  computed  as 

Ei/r  =  Ti  "  Hi  ‘  Tr  +  Hr  "  h/r  1  =  1»2,...,N 

where  all  quantities  have  previously  been  defineu.  The  maximum 
relative  error,  regardless  of  reference  station  bias,  from  this 
set  of  N  errors  is  then 

Emax  ^Ei/r^max  “  ^Ei/r  ^min 
at  each  grid  point. 

The  dashed-line  polygon  in  Figure  5  is  the  contour  of  the 
(known)  maximum  relative  error  (0.18  sec)  at  the  true  epicenter 
marked  at  X.  The  circle  enclosing  the  value  0.11  is  the  max¬ 
imum  relative  error  as  a  result  of  the  final  solution  obtained 
with  SHIFT-H61.  The  solid- line  polygon  is  the  contour  of  the 
maximum  relative  error  estimated  to  be  0.4  sec  (0.3  sec  higher 
than  that  at  the  final  solution)  for  this  event.  A  seismic 
analyst  can  usually  estimate  his  reading  errors  quite  well,  and 
if  the  effects  of  previously-determined  relative  travel-time 
anomalies  are  removed  (they  must  be  determined,  not  estimated), 
the  estimates  can  be  used  to  contour  the  maximum  relative  error. 
A  relative  error  estimate  of  0.4  sec  is  liberal  for  an  event  of 
the  size  of  Bourbon,  but  it  is  an  example  of  the  manner  of  using 
any  estimate.  The  approximate  area  of  the  estimated  polygon  is 
less  thar  340  km  .  The  area  of  the  corresponding  ellipse  pre¬ 
viously  computed  (Table  VII)  is  7118  km2. 
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Figure  6  shows  similar  results  for  Scotch.  Again,  an  error 
estimate  0.3  sec  higher  than  that  obtained  from  the  final  solu¬ 
tion  (0.5  sec)  is  used  to  contour  0.8  sec.  The  approximate  area 

2 

of  the  estimated  polygon  is  less  than  400  km  compared  to  2309 
km  of  the  standard  confidence  ellipse. 

Although  Scotch  was  a  large  event  and  the  reading  errors 
were  expected  to  be  small,  an  actual  maximum  relative  error  of 
0.7  sec  is  observed.  The  reason  for  this  appears  to  be  the 
slight  instability  of  the  relative  anomalies  for  the  Nevada  Test 
Site.  Scotch  was  located  about  40  km  to  the  northwest  of 
Bourbon,  and  hence  about  40  km  from  the  area  of  the  reference 
events  used  for  determining  the  anomalies.  For  networks  the 
size  of  that  used  for  Scotch  (80°  azimuth  aperture),  it  is  not 
surprising  that  the  anomalies  are  not  constant. 

Figures  7  and  8  show  the  maximum-relative-error  grid  for 
a  network  of  three  stations  obtained  by  deleting  SV3QB  from 
Bourbon  and  SV3QB  and  PG-BC  from  Scotch.  For  three  stations, 
it  is  impossible  to  estimate  the  location  errors  by  the  usual 
means  as  the  statistical  degrees  of  freedom  are  reduced  by  zero. 
However,  with  maximum  relative  errors,  reliable  estimates  can 
be  made.  Again,  for  Bourbon,  which  is  located  within  the  area 
of  the  reference  events,  an  error  estimate  of  0.4  sec  is  con¬ 
toured  (Figure  7),  and  for  Scotch,  40  km  from  the  area  of  the 
reference  events,  an  error  estimate  of  0.7  sec  is  used  (Figure 
8).  The  approximate  areas  of  the  polygons  for  Bourbon  and 
Scotch  are  350  km  and  630  km  respectively.  Therefore,  by 
using  estimates  of  the  maximum  relative  errors,  reasonable 
estimates  of  location  errors  can  be  made  when  the  usual  statis¬ 
tical  confidence  estimates  are  impossible  to  compute. 

Standard  Deviation  Ellipses. 

Subroutine  SIGRID  also  produces  an  output  of  network  stand¬ 
ard  deviations,  in  addition  to  the  maximum  relative  errors,  for 
a  similar  set  of  grid  positions. 
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Figure  9  shows  the  zero-mean  standard  deviation  output'  for 
Bourbon,  with  anomalies.  The  values  (X100)  at  each  grid  point  k 
are  computed  as 
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(Ak/r  ‘  Ai/r 


E  ) 
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where  A  is  the  computed  anomaly  at  k,  A  is  the  predetermined 
anomaly,  and  Er  is  the  average  error,  or  bias,  at  the  reference 
station  r.  This  output  actually  shows  how  the  least-squares 
procedure  within  SHIFT  minimizes  the  sum  of  squares  (or  standard 
deviation)  of  errors.  Tne  dashed-line  ellipse  is  the  contour  of 
the  known  standard  deviation  of  time  errors  for  Bourbon.  The 
standard  confidence  ellipse  would  be  the  same  as  that  shown  in 
Figure  5. 

Therefore,  by  using  the  parameters  of  the  computed  confi¬ 
dence  ellipse,  the  maximum-relative-error  SIGRID,  and  the 
standard-deviation  SIGRID,  the  utmost  information  is  being 
elicited  from  the  time  data. 
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IMPROVING  THE  ESTIMATES  OF  TRAVEL-TIME  ANOMALIES 

The  computed  travel-time  anomalies  used  in  this  report  were 
derived  from  a  few  selected  reference  events  to  show  the  tele- 
seismic  network  capability  for  locating  events  in  the  Nevada 
Test  Site  area.  It  is  known  that  the  anomalies,  for  some  stations 
and  for  this  area,  are  not  constant  but  exhibit  some  variability 
at  different  positions  within  the  area  shown  in  Figure  1.  A 
much  better  estimate  of  the  true  travel-time  anomaly  for  the  NTS 
area  can  be  obtained  by  averaging  the  anomalies  at  each  station 
for  all  19  events  in  this  study.  These  average  anomalies,  given 
in  Table  VUI  using  the  Herrin  66  travel-time  table,  should  be 
used  for  locating  new  events  in  the  NTS  area  when  a  network  is 
comprised  of  any  of  the  listed  stations.  Table  VIII  is  a  com¬ 
puter  output  of  program  TIMEANOM;  the  key  circled  numbers  shown 
are  described  in  Appendix  II.  Tables  IX  and  X  list  the  station 
anomalies  for  the  Jeffreys-Bullan  and  Herrin  61  travel-time 
tables  respectively. 

Other  Variables  for  Computing  Anomalies. 

All  of  the  station  anomalies  in  this  report  were,  of  course, 
computed  from  nuclear  explosions,  the  positions  of  which  are 
extremely  well  known  in  the  three  dimensions  of  latitude,  long¬ 
itude,  and  depth.  Past  studies  (Chiburis  and  Dean,  ibid,  and 
Chibuns ,  ibid)  and  studies  currently  in  progress  show  that 
earthquakes  are  not  quite  so  well  behaved  as  explosions  for  sev¬ 
eral  reasons.  First,  simple  epicenter  mislocations  can  yield 
errors  in  anomaly  estimates  as  high  as  two  or  three  seconds  for 
large  networks.  The  mislocation  effect  is  described  in  Chiburis 
and  Dean,  ibid,  p.  31  ff.  Second,  and  perhaps  more  serious, 
depth  effec:  a  must  be  taken  into  account.  There  a^e  no  reasons 
to  suppose  that  the  anomalies  measured  for  shallow  events  in  a 
particular  region  may  not  change  significantly  for  deep  events 
in  the  same  region.  Added  to  these  effects  are  the  problems  of 
depth  errors  in  the  located  events. 
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One  of  the  current  studies  at  the  SDL  has  shown  conclusively 
that  for  a  North  American  Network  of  21  i,RSM  stations,  the  station 
anomalies  computed  from  the  explosion  LONG  SHOT  on  Amchitka  Is¬ 
land  are  not  in  as  good  agreement  as  expected  with  the  anomalies 
computed  from  earthquakes  in  the  large  region  of  the  Rat-Andreanof 
Is.  at  depths  of  15  km  to  300  km  (discrepancies  as  large  as  1.5- 
2.0  sec).  At  least  part  of  the  difficulty  seems  to  be  associated 
with  hypocenter  mislooation,  because  when  the  set  of  earthquakes 
was  relocated  by  using  the  LONG  SHOT  anomalies,  most  of  the 
epicenters  shifted  an  average  of  about  50  km  after  which,  in 
general,  the  serious  time  errors  were  significantly  reduced. 
Certainly  the  LONG  SHOT  anomalies  are  not  valid  for  all  of  the 

events  in  the  set  because  the  linear  size  of  the  region  is  about 

800  km  or  more,  much  too  large  for  a  single  set  of  anomalies  to 
exist  throughout,  but  they  should  be  valid  for  events  very  near 
LONG  SHOT  and  these  still  shift  40-50  km.  Figure  10  shows  the 

directions  of  the  shifts,  indicating  that  large  bias  effects 

may  exist  for  the  Andreanof  locations  as  reported  by  the  U.S.C. 
and  G.S.  This  study  is  continuing. 

Anomalies  vs.  Residuals. 

The  preceding  results  indicate  that  comparable  location 
accuracies  for  explosions  are  obtained  when  using  either  rela¬ 
tive  anomalies  or  absolute  residuals.  The  question  may  arise 
then  as  to  what  differences,  if  any,  one  may  expect  when  trying 
to  determine  a  true,  physical  time  correction  for  an  earthquake 
region.  Some  of  the  possible  differences  are  as  follows: 

1.  Origin  time  errors  (serious  for  some  earthquakes)  play 
no  role  when  relative  times  are  used;  time  bias  is  removed,  so 
that  events  in  the  same  region  may  be  compared  with  one  another 
in  determining  the  actual  anomaly. 

2.  First  extrema  as  well  as  first  motions  can  be  read  to 
determine  consistent  corrections  from  different  events  in  the 
same  regions.  This  is  important  if  one  has  only  a  few  small 
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events  for  time  calibrating. 

3.  When  using  residuals,  the  defined  regions  can  be  greater 
in  number  and  smaller  in  area  than  when  using  anomalies,  and  the 
problem  of  time  calibrating  the  earth  becomes  more  difficult.  A 
simple  case  qualitatively  illustrates  this  point  where  Figure  11 
shows  a  geologically  faulted  area  composed  of  two  different  crus¬ 
tal  media  I  and  II  with  average  P-wave  velocities  of  V.^  and  V 2* 
The  network  stations  are  labeled  1,  2,  and  3  each  having  a  res¬ 
idual  from  an  event  in  medium  I  computes  an  R*,  R*,  and  R*  and 

anomalies  relative  to  station  1  computed  as  0,  A*^,  A*^.  For 
an  event  occurring  in  region  II,  however,  each  of  the  residuals 

have  a  bias  added  to  them  due  to  the  different  velocity  in 
the  vicinity  of  the  source  in  this  region;  the  anomalies,  how¬ 
ever,  would  be  expected  to  change  only  slightly  because  the 
hypocenter-station  raypath  differences  between  the  two  (tele- 
seismic)  sources  would  be  negligible.  Hence,  for  each  of  the 
stations  there  would  have  to  be  two  residual  regions  but  only 
one  anomaly  region.  Of  course,  for  stations  situated  on  the  side 
of  region  II  and  opposite  to  the  reference  station,  there  would 
then  be  two  regions  for  both  residuals  and  anomalies.  However, 
at  least  some  reduction  in  the  time-calibration  regionalization 
should  be  realized  when  relative  anomalies  are  used. 
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CONCLUSIONS 


The  following  conclusions  are  made  concerning  the  results 

obtained  by  analyzing  19  nuclear  explosions  detonated  within  a 
2 

2500  km  area  of  the  Nevada  Test  Site.  Seismograms  were  used 
from  teleseismic  stations  forming  networks  of  four  to  thirteen 
stations . 

1.  For  limited-station  teleseismic  networks,  the  location 
capability,  for  17  events  occurring  in  the  Nevada  Test  Site  and 
without  applying  previously-determined  residuals  or  travel-time 
anomalies,  is  about  26  km,  regardless  of  the  program  used  (LOCATE 
or  SHIFT)  and  regardless  of  the  travel-time  table  employed 
(Herrin  61,  Herrin  66,  or  JB  tables). 

2.  The  location  capability  for  the  same  networks  and  time 
data  is  better  than  3  km  when  previously-determined  residuals  or 
travel- time  anomalies  are  applied,  regardless  of  program  and 
travel- time  table.  That  is  to  say,  when  these  travel-time  curves 
can  equivalently  be  replaced  by  observed  travel  times  from  each 
station  to  an  accurately  known  explosion  point,  the  particular 
curve  selected  as  a  standard  is  essentially  irrelevant  and  epi¬ 
centers  of  other  nearby  explosions  can  be  located  within  the  in¬ 
dicated  error  limit. 

3.  The  effect  of  the  number  of  stations  in  the  range  from 
3  or  4  to  13  on  the  location  capability  of  the  networks  used  in 
the  study  is  negligible,  either  with  or  without  travel-time  anom¬ 
alies  . 

4.  The  effect  of  the  range  of  epicentral  distances  (distance 
aperture)  on  the  capability  of  locating  known  surface  events  with 
the  same  networks  is  negligible,  either  with  or  without  travel¬ 
time  anomalies. 

5.  The  azimuth  aperture  has  an  observable  effect  on  the 
location  capability  of  the  networks  where,  for  apertures  down  to 
60° ,  the  location  errors  without  travel-time  anomalies  are  as 
high  as  60  km,  and  with  anomalies  6  km. 
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6.  The  areas  of  computed  confidence  ellipses  can  be  reduced 
by  factors  of  1/5  to  1/152  with  the  application  of  travel-time 
anomalies  and  still  enclose  the  event  locations  which,  for  this 
study,  are  accurately  known,  as  are  the  explosions  used  to  cali¬ 
brate  the  praticular  epicenter-station  paths  used  in  the  study. 

In  effect,  these  reduced  ellipses  represent  uncertainties  in  the 
position  of  the  epicenters  relative  to  the  positions  of  the  calib¬ 
ration  events.  - 

7.  With  as  few  as  4  or  5  stations,  computed  confidence 
ellipses  with  anomalies  are  unrealistically  large  due  to  the  small 
number  of  statistical  degrees  of  freedom.  (The  Bourbon  explo¬ 
sion,  with  4  stations,  has  an  ellipse  7,000  km2  and  a  location 
error  of  1.9  km).  Estimates  of  the  maximum  relative  time  errors 
for  a  network  of  stations  permits  more  realistic  confidence  limits 
to  be  set  on  the  solution.  In  this  way,  Bourbon's  confidence 
polygon  is  about  340  km2.  Alsc,  as  few  as  three  stations  can  be 
used  to  obtain  the  confidence  limits  when  maximum  relative  errors 
are  estimated. 

8.  Travel-time  anomalies,  computed  from  a  few  selected 
events,  show  some  variability  across  the  NTS  area  in  Figure  1, 
but  the  effect  on  relative  location  accuracy  is  small.  The  size 
of  the  area  in  Figure  1  is  about  2500  km2,  implying  That  fairly 
large  regions  are  involved  for  determining  reasonably  constant 
anomalies,  so  the  problem  of  time  calibrating  a  stationary  net¬ 
work  for  any  region  of  the  earth  from  which  either  explosions 

at  accurately  known  locations  or  earthquakes  (bias  effects  not 
included)  are  recorded  is  not  formidable.  In  either  case,  if 
accurate  locations  are  not  independently  known,  the  epicentral 
solutions  with  anomalies  included  reduce  to  locations  relative 
to  one  another  with  the  actual  error  of  the  whole  set  remaining 
unknown . 
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9.  Relative  anomalies  for  earthquake  regions  can  be  simpler 
to  assess  than  absolute  residuals  because  (a)  origin  time  errors 
are  eliminated;  <b)  first-extremum  anomalies  and  first  motion 
anomalies  from  several  events  can  be  combined  (except  when  ob¬ 
vious  period  differences  are  noted);  and  (c)  depending  on  the 
region  and  on  the  network  geometry,  the  regions  for  which  the 
calibrations  need  to  be  determined  can  be  fewer  in  number  and 
larger  in  area. 
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TABLE  I 


Event  Information 


Event 

Date 

Origin  Time 

Latitude 

Longitude 

No.  of 
Stations 

AUK 

02 

Oct 

64 

20 

03 

00.0 

37.078N 

116.009W 

6 

FORE 

16 

Jan 

64 

16 

00 

00.1 

.142 

.049 

12 

DUMONT 

19 

May 

66 

13 

56 

28.1 

.111 

.058 

8 

CHARTREUSE 

06 

May 

66 

15 

00 

00.1 

.348 

.322 

7 

TURF 

24 

Apr 

64 

20 

10 

00.2 

.150 

.055 

10 

KLICKITAT 

20 

Feb 

64 

15 

30 

00.1 

.151 

.040 

9 

PILEDRIVER 

02 

Jun 

66 

IS 

30 

00.1 

.  227 

.055 

7 

BRONZE 

23 

Jul 

65 

17 

00 

00.0 

.098 

.033 

9 

CORDUROY 

03 

Dec 

65 

15 

13 

02.1 

.165 

.052 

9 

BUFF 

16 

Dec 

65 

19 

15 

00.0 

.073 

.029 

8 

GREELEY 

20 

Dec 

66 

15 

30 

00.1 

.302 

.408 

13 

PIRANHA 

13  May 

66 

13 

30 

00.0 

.087 

.034 

9 

NASH 

19 

Jan 

67 

16 

45 

00.1 

.144 

.136 

7 

BOURBON 

20 

Jan 

67 

17 

40 

04.1 

.100 

.004 

4 

AGILE 

23 

Feb 

67 

18 

50 

00.0 

.127 

.067 

7 

COMMODORE 

20 

May 

67 

15 

00 

00.2 

.130 

.064 

5 

SCOTCH 

23 

May 

67 

14 

00 

00. 0 

.275 

.  370 

5 

BILBY 

13 

Sep 

63 

17 

00 

00.1 

.061 

.022 

14 

TAN 

03 

Jan 

66 

14 

00 

00.0 

.068 

.035 

8 

TABLE  II 


Station  Anomalies  and  Residuals 


Anomaly  Relative  to 

RK-ON 

Residual 

Event 

Station 

Herrin  61 

Herrin 

66  JB 

Herrin  61 

Key* 

AD- IS 

+1.11  sec 

+0.04  sec  +0.57  sec 

-1.1  sec 

B 

AX2AL 

+  2.21 

+  1.61 

+  2.23 

+  0.2 

A 

BE-FL 

+  1.11 

+  1.11 

+  1.15 

-0.9 

A 

BL-WV 

+  0.66 

+  0.61 

+  0.84 

-1.6 

A 

BR-PA 

+  0.74 

+  0.76 

+  0.84 

-1.5 

A 

CP008 

+  1.51 

+  0.86 

+  1.55 

-0.6 

A 

DH-NY 

+  0.09 

+  0.03 

+  0.11 

-2.2 

A 

EB-MT  ' 

-0.84 

+  0.08 

-0.68 

-3.1 

A 

EN-MO 

+  0.08 

+  0.54 

+  0.04 

-2.0 

B 

EU-AL 

+  2.66 

+  1.94 

+  2.71 

+  0.6 

A 

GG-GR 

+  1.06 

+  1.00 

+  2.07 

-1.2 

A 

HN-ME 

+  1.11 

+  0.55 

+  0.97 

-1.0 

A 

KC-MO 

+  1.71 

+  1.81 

+  1.45 

-0.7 

D 

LV-LA 

+  1.51 

+  1. 90 

+  1.48 

-0.6 

A 

LZ-BV 

+  0.90 

+  0.94 

+  1.19 

-1.9 

A 

NP-NT 

+  1.88 

+  0.96 

+  1.52 

-0.2 

A 

00-NW 

-0.01 

+  0.01 

+  0.52 

-2.4 

A 

PG-BC 

+  2.66 

+  2.94 

+  2.57 

+  0.3 

C 

PZ-PR 

+  2.28 

+  1.17 

+  1.  53 

+  0.1 

A 

SI-BC 

+  1.87 

+  2.81 

+  2.  05 

-0.2 

A 

SV2QB 

+  0.35 

-0.37 

+  0.16 

-2.0 

A** 

SV3QB 

+  0.35 

-0.37 

+  0.16 

-2.0 

A 

WH2YK 

+  1.36 

+  1.13 

+  1.50 

-0.8 

C 

RK-ON 

0 

0 

0 

-2.1 

A 

*A  =  Bilby-Tan 
B  =  Bronze- Corduroy- Buff 
C  =  Nash-Agile-Commodore 
D  =  Gree  ey 


**Set  equal  to  anomaly  for  SV3QB  determined  from  A 


TABLE  III  (Continued) 


TABLE  IV 


Location  Errors  when  Neither  Residuals 


nor  Anomalies 

are  Used 

(EVENT 

NO.  OF 
STATIONS 

LOCATE -H61 

SHTFT-H61 

SHIFT-H66 

SHIFT-JB 

Buff 

8 

9.1 

8.6 

10.9 

5.7 

Turf 

10 

14.  5 

15.0 

8.0 

18.3 

Corduroy 

9  . 

11.6 

10.8 

5.2 

6.3 

Nash 

7 

17.2 

18.1 

13.4 

23.1 

Bourbon 

4 

4S.0 

49.0 

17.6 

42.6 

Piledriver 

7 

39.6 

39.3 

37.7 

41.9 

Bronze 

9  < 

13.3 

13.1 

8.5 

12.3 

Auk 

6 

7.1 

7.2 

3.7 

4.8 

Piranha 

9 

41.1 

40.4 

45.5 

46.3 

Fore 

12 

38.1 

38.5 

20.0 

36.6 

Greeley 

13 

10.1 

12.3 

11.6 

25.5 

Dumont 

8 

43.0 

43.6 

44.9 

47.1 

Chartreuse 

7 

60.0 

59.8 

48.5 

62.8 

Commodore 

5 

11.1 

10.9 

23.3 

12.6 

Agile 

7 

23.4 

23.0 

15.0 

27.7 

Scotch 

5 

17.4  * 

17.8 

26.6 

19.6 

Klickitat 

9 

33.5 

33.7 

12.8 

38.2 

l 

438.1 

441.1 

353.2 

471.4 

Mean  error ,  km 

25.8 

25.9 

20.8 

27.7 

tw*we» 


TABLE  V 


Location  Errors  When  Residuals 
or  Anomalies  are  Used 


EVENT 

STATIONS 

LOCATE -H61 

SHIFT-H6I 

SHIFT  H6 6 

SHIFT 

Buff 

8 

2.4 

0.1 

0.3 

0.2 

Turf 

10 

0.5 

0.1 

0.2 

0.1 

Corduroy 

9 

1.6 

1.0 

0.9 

0.9 

Nash 

7 

5.8 

1.7 

1.4 

1.8 

Bourbon 

4 

2.8 

1.9 

2.1 

1.9 

Piledriver 

7 

2.0 

2.9 

1.5 

3.3 

Bronze 

9 

4.2 

3.0 

2.7 

2.9 

Auk 

6 

3.1 

3.0 

2.9 

3.0 

Piranha 

9 

0.8 

3.1 

2.7 

3.2 

Fore 

12 

3.3 

3.2 

2.9 

3.1 

Greeley 

13 

1.0 

3.3 

3.8 

3.3 

Dumont 

8 

2.2 

3.4 

3.0 

3.4 

Chartreuse 

7 

5.  7 

3.6 

3.7 

4.2 

Commodore 

5 

4.1 

3.7 

3.2 

3.3 

Agile 

7 

1.5 

3.9 

2.4 

3.1 

Scotch 

5 

3.0 

5.0 

4.8 

5.7 

Klickitat 

1 

9 

6.1 

6.0 

5.6 

6.2 

50.7 

48.9 

44.1 

49.6 

Mean  error, 

km 

2.98 
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2.59 
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TABLE  VI  (Continued) 


TABLE  VII 


Confidence  Regions 
SHIFT-61 


Event 

Ellipse 

Without  anomalies 
a 

Areas  km 

With  anomalies 
b 

Improvement 
Factor,  a/b 

Nc  of 

S  t n  t i ons 

Fore 

2059* 

348 

6 

12 

Buff 

2402 

50 

48 

8 

Chartreuse 

1041* 

131 

8 

7 

Auk 

10035 

66 

152 

6 

Bourbon 

206086 

7118 

29 

4 

Dumont 

1495 

52 

29 

8 

Agile 

4020 

234 

17 

7 

Nash 

7256 

113 

64 

7 

Commodore 

27698 

288 

96 

5 

Greeley 

2620 

96 

27 

13 

Klickitat 

2384 

202 

12 

9 

Turf 

2513 

524 

5 

10 

Piranha 

2159 

15 

144 

9 

Scotch 

37869 

2309 

16 

5 

Corduroy 

15S4 

25 

62 

9 

Bronze 

2038 

25 

81 

9 

Piledriver 

1593 

79 

20 

7 

^Confidence  ellipse  does  not  contain  true  epicenter 
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TAliir  VIII 

«  N  A  V  E  l  •  T  |  N  E 


A  N  0  N  A  l  I  ►  5 
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Figure  1.  Nevada  Test  Site  Area 
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Figure  2.  Location  error  vs  number  of  recording  stations 
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Figure  3.  Location  error  vs  azimuth  aperture 
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Figure  5.  Maximum-relative-error  output  from  SIGRID;  event  Bourbon 
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Figure  10.  Location  Shifts  with  LONGSHOT  Anomalis 


Figure  11.  Faulted  region  illustrating  difference 
when  time-calibrating  with  residual  R 
or  anomolies  A. 


APPENDIX  I 

INPUT  SEQUENCE  FOR  PROGRAM  SHIFT 


Figure  A-l 


INPUT  CARD  FORMATS  AND  PARAMETERS 
OPTION  CARD  -  FORMAT  (7I2.4FS.0,  413) 

NATS  -  Number  of  stations  and  anomalies  in  entire  run  of 
several  cases.  (Maximum  of  50) 

=  0:  set  to  21  internally 

LSW  =  1:  Read  input  data  beginning  with  travel-time 
table  (See  Figure  A), 

=  2:  Read  input  data  beginning  with  station  deck. 

=  3:  Read  input  data  beginning  with  anomaly-heading 
card. 

=  0  or  4:  Read  input  data  beginning  with  epicenter 
card. 

NJUMP  =  0:  Read  J-B  table  (Tape  940,  SDL) 

=  1:  Read  Herrin  66  table  (Tape  940,  SDL). 

LSC  =  0:  Read  anomalies  in  Format  A,  arrival  times  in 
Format  B. 

=  1:  Read  anomalies  in  Format  A,  arrival  times  in 
Format  C. 

=  2:  Read  anomalies  in  Format  D,  arrival  times  in 
Format  E. 

IJK  =0:  No  SIGRID  computation 

i  0:  Compute  SIGRID 

ISKIP  =  0:  Program  locates  epicenter 

i  0:  Program  does  not  locate  epicenter;  computes 
SIGRID  (if  IJK  t  0)  using  input  location 
coordinates. 

NOUT  =  Flag  requesting  maximum-relative-error  station  sub¬ 
scripts  in  SIGRID 

=0:  Do  not  print  subscript  matrix. 


-Al- 


t  0:  Print  subscript  matrix  . 

CT  =  Distance  increment  for  computing  time  derivatives 

=  0:  Set  to  0.01°  arc  internally. 

SLOW  =  Factor  to  reduce  the  size  of  the  iterative  corrections 
to  the  event  parameters 

=  0:  Set  to  1  internally. 

TEST  -  Incremental  distance  criterion  for  convergence 

=  0:  Set  to  0.0001  internally. 

CLAT  =  Grid  increment,  kilometers,  for  SIGRID  computation 

=  0:  Set  to  1.0  internally. 

NCOL  =  Number  of  columns  from  SIGRID  output  grid 

=  0:  Set  to  30  (maximum)  internally, 

KKEND  =  Number  of  lines  from  SIGRID  output  grid 

s  0:  Set  to  30  (maximum  of  100)  internally. 

ITER  s  Number  of  iterations  allowed  for  location  con¬ 
vergence 


=  0:  Set  to  10  internally. 

IMK  =  0:  Use  input  arrival  times  in  SIGRID 

t  0:  Use  predicted  arrival  times  to  grid  center 
in  SIGRID  . 


STATION  CARDS  - 


FORMAT  (A5 ,  3X,  2(FS.O,  F3.0, 


F5.1,  Al),  24X ,  F5.0) 


Column 


Parameter 


1-5 

9-13 

14-16 

17-21 

22 

23-27 

28-30 

31-35 

36 

61-65 


Station 

Station 

Station 

Station 

N  or  S 

Station 

Station 

Station 

E  or  W 

Station 


name 

latitude,  degrees 
latitude,  minutes 
latitude,  seconds 

longitude,  degrees 
longitude,  minutes 
longitude,  seconds 


elevation,  meters  (not 


necessary) 


X 
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ANOMALY -HEADING  CARD  -  FORMAT  (10A8) 


Description  of  anomalies,  region  name,  etc.  May  be  blank 


card . 


ANOMALY  CARDS 


Format  A  (8F10.4) 

Format  D  (10X,  F10.3) 

Anomalies  must  be  punched  in  same  order  as  the  station  cards 
and  equal  to  the  number  of  stations .  May  be  blank  cards  if  un¬ 
known  . 

EPICENTER  CARD  -  Format  (2A8,  A4,  2(F9.3,  Al),  F5.0) 


Column 

1-20 

21-29 

30 

31-39 

40 

41-45 


Parameter 

Arbitrary  name  of  event 
Estimated  latitude,  decimal  degrees 
N  or  S 

Estimated  longitude,  decimal  degrees 
E  or  W 

Estimated  depth,  kilometers 


ARRIVAL-TIME  CARDS  -  (Zero  input  arrival  time  indicates  no  read¬ 
ing  for  that  station) 

FORMAT  B  ( 8 ( 2F2 . 0 ,F6 . 3 ) ) 


Column 

1-2 

3-4 

5-10 

11-12 

13-14 

15-20 


Parameter 

Arrival  hour,  station  1 
Arrival  minute,  station  1 
Arrival  second,  station  1 
Arrival  hour,  station  2 
Arrival  minute,  station  2 
Arrival  second,  station  2 


Repeat,  eight/card,  until  no.  of  times  =  no.  of  stations 
FORMAT  C  (7F10.4) 

Fields  of  10,  seven/card.  Arrival  time  in  seconds  only. 
FORMAT  E  (A5 ,  5X,  2F2.0,  F6.3) 
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If  this  format  is  used,  one  card,  Format  (15),  indicating 
number  of  arrival  times  to  be  input,  is  required  prior  to  time 
cards  which  follows : 


Column 

Parameter 

I- 5 

II- 12 
13-14 
15-20 

Station 

Arrival 

Arrival 

Arrival 

name 
time , 
time , 
time , 

hour 

minute 

second 

With  this  format,  one  card/station  with  an  arrival  time  is  re¬ 
quired. 
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APPENDIX  II 


Description  of  TIMEANOM  Output. 

The  following  description  explains  the  presentation  of  the 
results  in  the  computer  output  of  the  TIMEANOM  with  these  refer¬ 
ence  numbers  appearing  in  Table  VIII. 

1.  Source  of  expected  travel  times.  For  Table  VIII,  the 
Herrin  table,  November  1966  version,  is  used;  for  Table  IX,  the 
JB  table;  for  Table  X,  the  Herrin  61  table. 

2.  Reference  station,  R,  selected  for  computing  relative 
anomalies.  In  this  report,  all  anomalies  are  relative  to  station 
RK-ON.  The  following  relation  may  be  used  to  change  reference 
stations ; 


A •  /  •  —  A .  ,  —  A.  , 

1/3  l/r  ]/r 

where  A^j  is  the  anomaly  at  station  i  relative  to  a  new  reference 
station  j . 

3.  All  expected  travel-times  in  this  report  have  been  cor¬ 
rected  for  the  ellipticity  of  the  earth  such  that  the  computed 
anomalies  may  be  used  in  conjunction  with  other  programs  requiring 
these  corrections. 

4.  An  arbitrary  geographic  name  given  to  the  event  region. 

5.  Range  of  epicentral  distance  in  the  event  region. 

6.  Range  of  epicentral  azimuth  in  the  event  region. 

7.  Date  and  arbitrary  name  given  to  each  event. 

8.  Epicentral  distance,  in  kilometers,  from  the  reference 
station,  R. 

9.  Epicentral  azimuth,  in  degrees  measured  from  north  to 
east,  from  the  reference  station,  R. 

10.  Station  designator,  i. 
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11.  Measured  travel- time  anomaly  in  seconds,  at  station  i 
relative  to  station  R  for  the  kth  event; 


Hk  +  Hk 
1  r 


where  T  is  the  observed  arrival  time  and  H  is  the  expected  (Herrin 
1966)  travel  time  from  the  hypocenter  of  the  kth  event  including 
correction  for  ellipticities  but  not  for  station  elevations. 

12.  A  fixed-point  zero  anomaly  indicates  that  no  reading 
was  made  at  the  station  for  that  event. 


13.  The  average  anomaly  at  station  i  of  N  recorded  events; 

N 


A., 

l/r 


l  *1 


,k=l 


i/r 


/N 


for  the  defined  region. 

14.  Standard  deviation,  or  error  of  estimate,  at  the  ith 
station  for  N  observations: 

1/2 


IT  N 

•v 

tu  “k-w* 

/(N-l)> 

for  the  defined  region. 

15.  Number  of  observations,  N,  at  station  i  for  the  defined 
region. 

16.  Total  number  of  epicenters  included  in  the  defined  re¬ 
gion. 

17.  Epicenter  latitude,  degrees  (USC6GS ) ;  plus  north,  minus 
south. 


18.  Epicenter  longitude,  degrees  (USC6GS) ;  plus  east;  minus 
west. 


19.  Event  depth,  kilometers  (USC6GS) . 

20.  Event  origin  time,  hours,  minutes,  seconds  (USC6GS). 

21.  Standard  deviation,  or  error  of  estimate,  of  the  kth 
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event  in  the  defined  region; 


o 


k 


1/2 


where  L  is  the  number  of  stations  recording  the  kth  event  not  in¬ 
cluding  the  reference  station  R. 

22.  Average  error,  or  bias,  of  the  kth  event; 

Ek  =i\1  (Ai/r  '  Si/r)/L 

where  L  is  the  number  of  stations  recording  the  kth  event  not 
including  the  reference  station  R. 

23.  Standard  deviation  of  the  kth  event  in  the  defined 
region,  with  the  reference-station  bias  removed: 


24.  Number  of  stations,  L,  recording  the  kth  event,  not 
including  the  reference  station  R. 

The  program  TIMEANOM  permits  a  rapid  determination  of  travel¬ 
time  anomalies  for  a  network  and  for  a  set  of  events  within  a  re¬ 
gion,  and  it  can  be  used  to  isolate  spurious  readings  or  blunders 
at  the  stations  or,  for  earthquakes,  possible  mislocations . 
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